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Abstract. The Yablonskii-Vorob'ev polynomials y n (i), which are defined by a second 
order bilinear differential-difference equation, provide rational solutions of the Toda 
lattice. They are also polynomial tau-functions for the rational solutions of the second 
Painleve equation (P77). Here we define two- variable polynomials Y n (t, h) on a lattice 
with spacing h, by considering rational solutions of the discrete time Toda lattice as 
introduced by Suris. These polynomials are shown to have many properties that are 
analogous to those of the Yablonskii-Vorob'ev polynomials, to which they reduce when 
h = 0. They also provide rational solutions for a particular discretisation of P77, 
namely the so called alternate discrete P77, and this connection leads to an expression 
in terms of the Umemura polynomials for the third Painleve equation (P777). It is 
shown that Backlund transformation for the alternate discrete Painleve equation is a 
symplectic map, and the shift in time is also symplectic. Finally we present a Lax pair 
for the alternate discrete P77, which recovers Jimbo and Miwa's Lax pair for P77 in 
the continuum limit h — ► 0. 
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1. Introduction 



The Toda lattice 

d 2 x, 



n 



n E Z 



(1) 



dt 2 



— e 



was the first integrable differential-difference equation to be discovered [Tj. The 
Yablonskii-Vorob'ev polynomials (2J [3] yield rational solutions of both the Toda lattice 
and the second Painleve transcendent (Pn), since the tau-functions of Pjj satisfy the 
bilinear form of the Toda lattice. In a previous work [I] one of the authors obtained an 
expression for solutions of the Toda lattice as ratios of Hankel determinants, by using 
the associated Lax pair to construct continued fraction solutions to a sequence of Riccati 
equations. This in turn led to an expression for the Yablonskii-Vorob'ev poynomials as 
Hankel determinants [5], equivalent to that discovered more recently j6] (see also [7]). 

Here we will consider the case when the time evolution becomes discrete. Our 
approach is to start from the Lax pair for the discrete Toda lattice (dTL) given by Suris 
[8]. In section 2 we present this construction and derive polynomials Y n (t,h) in two 
variables t, h, which satisfy relations at a discrete set of points t = t + mh, m e Z 
(where to is arbitrary). These polynomials tend to the Yablonskii-Vorob'ev polynomials 
y n (t) as the spacing h — > 0, so that Y n (t,0) = y n (t). A discrete analogue of the bilinear 
defining equation for the y n is derived and a corresponding representation of these Y n 
as Hankel determinants is given. 

The first few Yablonskii-Vorob'ev polynomials are 



A further property of these polynomials, investigated recently by Clarkson and Mansfield 
[9], is the distribution of their zeroes. It is known that each y n (t) has no zeroes in 
common with y n ±i(t), and that these zeroes are simple pm E]- Numerical studies 
indicate that these zeroes lie in approximately triangular arrays and that the zeroes of 
y n (t) interlace, in a certain sense, with those of y n+ i(t), in a similar way to the zeroes 
of classical orthogonal polynomials. In section 3 we show that the polynomials Y n (t, h) 
have analogous properties. 

In the continuum case the bilinear differential-difference equation that defines the 
Yablonskii-Vorob'ev polynomials is given in terms of the Hirota derivative D t by 



with the initial polynomials yo = 1, y\ = t, and this equation follows from the 
representation of Pjj as a pair of bilinear equations for the associated tau-functions. In 
section 4 we consider a discretisation of Pu, namely the alternate discrete Pn equation 
(alt-dPjj) studied in [12], and explain how its solutions are specified by tau-functions 
which satisfy a discrete bilinear equation together with a quadrilinear (degree four) 
relation. The bilinear relation defining the polynomials Y n (t, h) is a consequence of 
these tau-function equations. The alt-dPji equation has Backlund transformations 
corresponding to the shifts n — > nil, just as in the continuum case, which imply 



yo = 1, yx = t, y 2 = t 3 + 4, y 3 = t e + 20t 3 - 80, y 4 = t 10 + 60t 7 + 11200*. 



(2) 



yn+lVn-l = tyl - 2£>i Vn ' V; 



(3) 
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that the discrete Yablonskii-Vorob'ev polynomials satisfy recurrence relations which 
are analogous to those in the continuum setting. It is also known that the alt-dPjj 
equation arises as the contiguity relations for a Backhand transformation for Pin- The 
latter connection leads to an alternative formula for the discrete Yablonskii-Vorob'ev 
polynomials in terms of determinants of Jacobi-Trudi type, corresponding to Umemura 
polynomials [IB"] . 

It is known from the work of Okamoto that the continuum Pjj can be written 
in Hamiltonian form, and the Backlund transformation is a canonical transformation 
[14] . In section 5 we present a Poisson structure for the discrete case such that 
both development in time and the Backlund transformation for alt-dPji are symplectic 
maps. Finally in section 6 we present a Lax pair for alt-dPjj, which relates it to the 
isomonodromic deformation of an associated linear system, and show that this tends to 
Jimbo and Miwa's Lax pair for the continuum case as h — > 0. 



2. Discrete time Toda lattice and discrete Yablonskii-Vorob'ev polynomials 



l+h 2 cxp(x„^ 1 (t)~x n (t)) 
1+h 2 cxp(x„(t)-x n+1 (t)) 



(4) 



The equations for the discrete Toda lattice obtained by Suris [8] are 
x n (t + h) - x n (t) = hn n (t + h), 
ir n (t + h) - 7T n (t) = ±log 

where 7r n denotes the canonically conjugate momentum to x n , and these discrete 
equations clearly yield Hamilton's equations for the Toda lattice (CEJ) in the continuum 
limit, as h — *■ 0. They arise from the consistency condition for the Lax pair 

tt n (i + h)= V n (t)V n (t), * n+ i(t) = L n (t)V n (t), (5) 

that is L„(t + h)"V n (t) = V n+ i(t)L„(t), where 

(f s ( Aexp(-/wr n (*)) - \ hexp(-x n (t)) \ ( ) 

A> \-hexp(x n (t)-hn n (t)) j {0) 

and 

V(t\-( i -hexp{-x n {t))\ 

Upon setting * n (t) = [X n (t),Y n (t)] T and Z n (t) = X n (t)Y n {t)~ l , we find from the 
second of (jS]) that 

z = hexp(-x n (t)) 

Aexp(-/i7r n (t)) - \ + hexp(-hiT n (t) + x n (t))Z n+1 (t) 
This recurrence relation can be used to generate the continued fraction expansion 

XhZ {t) = -Jl^ -ll^ -Jh— . , (9) 

where 7 = A 2 , fi = h 2 exp(— xo(t)), g\ = — exp(— hn (t)) and for n = 2, 3, . . . we have 
f n = ~ /i 2 exp(-x„_i(t) + x n - 2 (t) - hTr n - 2 {t)), 

g n = - exp(-/i7r ri _i(t)). (10) 
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Similarly from the first of §5$) we have the discrete time Riccati equation 

Z (t + h) - Z (t) = c -oW + " l) Zo(t) - je*-^Z (t + h)Z (t). (11) 

In our previous work [15] we took x_i{t) — > — oo so that this became a linear equation. 
Here we take x_i(i) = and set Z (t) = TrW(t) so that with 7 = A 2 we find 



Wit + h) = -h 2 e~ x °® + -Wit) - -Wit + h)W(t). 

7 7 



(12) 



Now the continued fraction ([9]) has expansions both in positive and negative powers 
of 7, namely Wit) = Y^=i a nit)^/ n and Wit) = Yl^Zo a nif)l n respectively. Substituting 
these expansions in turn into ( jl~2l) and equating coefficients of corresponding powers of 
7, we obtain the recurrence relations for their coefficients, given respectively by 

3 

= aj{t + h) + a k {t + h)aj^ k+1 (t), j = 1,2,3,..., (13) 



k=l 




-2 



(14) 



aj (t + h) = a j+1 (t) - 22 ® k (t + h)a j+1 - k (t), j = - 

k=j+i 

where ai(t) = h 2 exp(— xo(t)) and a^it + h) = — h 2 exp(— xo(t)). 

The continued fraction in is known as a T-fraction [16] and its elements are 
given in terms of Hankel determinants. If we take the definitions 

Ht n+2) , v n = Hi n+1 \ (15) 



u. 



with Hankel determinants 



H, 



(m) 



/3m 
An+1 



An+1 
An+2 



fim+k-l fim+k ■ 

whose elements are given by 

Pk = -a k (t), k = 1,2, 

= a k (t), k = 0,-1,-2,..., 
then the elements of the T-fraction are given by 

V n -2U n Vn-lUn 



m+2k-2 



(16) 



fn 



9n 



n 



2,3,4, 



(17) 



(18) 



From the expression for the f n ,gn in terms of the coordinates and momenta given 
by (TTU]) and the equation of motion (jlj), it follows that 



fi(t + h) 
fi(t) 

and more generally 
f n+1 (t + h) g n+1 (t + h) 



-9i{t + h), 



9i(t + h) 
9i(t) l + 



1 _i_ MQ. 

^ fli(t) 



/n+l(0 



0n(*) 



g n +i(t + /Q 

0n+l(*) 



1 + 



1 + 



U+2(t) 
9n + l(t) 

9n(t) 



tiGff. 



(19) 



(20) 
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Proposition 2.1 The two types of Hankel determinant defined by < T75]) are related by 
v n (t + h)=u n (t), n = 0,l,2,.... (21) 

Proof: Substituting the expressions for /„,, g n given by ( jTSl) into the first of ( 120]) one 
finds that 

^ n +i(t + ft> w (t) _ vi(t + h)u (t) 
v n (t + h)u n+ x(t) v (t + ft)wi(t)' 

But m = /Jq 2 - 1 = 1 = fTp = t> and Wi(t) = fix = — ft 2 exp(— x (£)), = /3 = 

— /i 2 exp(— x (t — h)). Therefore the right hand side of (12"2"]) equals unity and the result 
follows. □ 

Proposition 2.2 When exp(—xo(t)) = —4, i/ie Hankel determinants u n satisfy the 
bilinear relation 

h 2 t\ 

1 ^- J + ft)«n(£ - ft) = u„(t) 2 + u n+ i(t)u n _i(t) (23) 

/or n = 1, 2, . . .. 

Proof: Substituting for g n and / n from ([T5]) into the second of f[2"U]) and then eliminating 
the dependence on the v n using (I2~T|) . the result follows from the fact that Ux(t) = 
Vl (t + h) = ¥± and u 2 {t) = -g(t 3 - h 2 t + 4). □ 

Remark. In [T7j, Hankel determinant solutions are given for a different (but gauge 
equivalent) bilinear form of the discrete Toda lattice equation, namely the equation 
P^Pn 1 ~ (Pn) 2 = £ 2 Pn+iPn-v Upon setting s 2 = 1 and t = (I - n)h, u n {t) = P l +ht / P % 
satisfies (123]) for a suitable choice of the initial condition Pq* for n = 0. 

It follows from the recurrence relations (11411 that when exp(— xo(t)) — — 4 the ctj(t), 
and hence also u n (t) and v n (t), are polynomials in ft,, t. We now renormalise the u n so 
that they are O(h ) as h — > 0. 

Definition 2.1 T/ie discrete Yablonskii-Vorob'ev polynomials Y n (t, h) are defined by the 
bilinear recurrence 

h 2 Y n+1 (t, h)Y n _x(t, h) = (h 2 t - 4) y n (t + h, h)Y n (t -h,h) + AY n (t, h) 2 (24) 
for n > 1, with Yo(t, h) = 1, Yx(t, h) = t as initial data. 

Theorem 2.1 The discrete Yablonskii-Vorob'ev polynomials Y n (t,h) G Z[t, h 2 ] 
are given in terms of the Hankel determinants u n in ( f73j) by Y n (t,h) = 
{—l) n {— r -j-)~ n ( n+1 ^ 2 u n {t) for n = 0, 1, 2, . . .. Each Y n is a monic polynomial of degree 
n(n + l)/2 in t, satisfying Y n (t, h) = y n (t) + 0{h 2 ) as h — » 0, where the y n [t) are the 
usual Yablonskii-Vorob'ev polynomials defined by {^) with y — 1, yx — t. 

Proof: With ax(t) = —h 2 t/4: it is easy to prove by induction from the first recursion 
in flUJ) that a^t) = h 2 Pj/A j for j > 1, where Pj E Z[t,h]. Similarly, with a (t) = 
h 2 (t — h)/4 the second recursion in (JT3J) implies that otj(t) = h 2 Pj/ for j < 0, 
where Pj e Z[t, ft]. It follows from their definition in terms of the matrix elements (1171) 
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that the Hankel determinants u n are polynomials in Q[t, h], with powers of 4 as the only 
possible denominators of the coefficients. If we let Y n (t,h) = ( y —l) n (—^-)~ n( - n+1 ^ 2 u n (t) 
for iieN then we find Yo(t, h) = 1, Yi(t, h) = t, and then the relation ( 1241) follows from 
Proposition 12. 2\ by substituting for u n (t) in ff23l) in terms of Y n (t, h). Since the Y n are 
uniquely defined by the recurrence ( 124"1) together with the given initial data, the formula 
in terms of renormalised Hankel determinants guarantees that they are polynomials in 
t, and the recurrence also implies that they are monic and of degree d n := n(n + l)/2) 
in t. However, further analysis is required to verify that there are no powers of h or 
powers of 4 in the denominator for this choice of normalisation. 

By induction, suppose that Y n (t,h) G Z[t,h], and hence is regular as h — > 0, for 
n = 0, 1, . . . , N (which clearly holds for N — 1). Taking a Taylor expansion in t, with 
derivatives denoted by Kjv,(j)t = ^r"; we see that 

Y N ( t +h,h)Y N (t-h,h) = (e£ /2] fe2fcy (gr (u) ) 2 - h2 (sEr 1)/21 ^5^ ) 2 

= Y N (t,h) 2 + h 2 P{t,h), 

where P(t, /i) is a polynomial in t and /t by the inductive hypothesis. Moreover, at 
leading order we have 

P(t, h) = Y N (t, 0)Y Nitt (t, 0) - Y N , t (t, 0) 2 + 0(h). 

Substituting this into the right hand side of (|24l) and dividing by h 2 we have 

Y N+1 (t, /i)yjv_i(t, /i) = tY N (t + h, h)Y N (t -h,h)- 4P(t, /i), (25) 

from which it follows that Y^ + i(t, h) is regular as h — > 0, and hence lies in Q[t, /i] with 
at worst powers of 4 as denominators of its coefficients. Now for some K > we can 
write Y N+ i(t,h) = Y N+1 (t, h)/A K where Fjv+i £ %[t,h] with 4 /fF A r +1 . If we multiply 
both sides of (124]) through by 4^ we see that if K > then we have 4|/i 2 l7v + iY/v-i; but 
then since Yat-i is monic in t it is clear that 2 / Y/v-i(£, ^) we must have 4|Y/v + i, which 
is a contradiction. Hence K = and Y/v+i G Z[t, /i] as required. 

Setting /i — > — /i in the recurrence (124"]) it also follows immediately by induction 
that Y n (t, h) = Y n (t, —h) for all n, so in fact we have polynomials in Z[£, /i 2 ]. Thus we 
can write Y n (t, h) = y n (t) + 0(h 2 ) where y (t) = Y Q (t, h) = 1 and y x {t) = Yx(t, h) = t, 
and using the leading order part of P(t, h) in (1251) we find that for n > 1 the y n satisfy 

Vn+xVn-i = ty 2 n - Ay n y n + 4y 2 , 

which is precisely the defining relation ([3]) for the usual Yablonskii-Vorob'ev 
polynomials. □ 

The first few discrete Yablonskii-Vorob'ev polynomials are 

Y (t,h) = l, 
Y 1 (t,h)=t, 
Y 2 (t,h) =t 3 + 4-h 2 t, 

Y 3 (t, h)=t 6 + 20t 3 - 80 + h 2 (4t - 5t 4 ) + 4t 2 h\ 
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Y 4 (t, h) = t 10 + 60t 7 + 11200t - h 2 (15t 8 + 252t 5 + 3360t 2 ) + 

+ /i 4 (63t 6 + 480t 3 + 576) - /i 6 (85t 4 + 288t) + 36h 8 t 2 . 



(26) 



These examples clearly reduce to the usual Yablonskii-Vorob'ev polynomials ([2]) when 



3. Zeroes of discrete Yablonskii-Vorob'ev polynomials 

In the continuum case it has been shown that the zeroes of y n (t) are simple and are 
not zeroes of y n+ i(t) [101 [TTJ. Here we prove an analogous result in the discrete case, 
by considering the roots of Y n (t,h) as a polynomial in t, that is to = to(h) such that 
Y n (t (h), h) = 0. To begin with we require a simple observation. 

Lemma 3.1 The discrete Yablonskii-Vorob'ev polynomials never vanish at t = 4/h 2 . 
More precisely, 



Proof: The result is trivially true when n = 0, 1. Assume that it holds for n = 
1,2, ...,N. Then from ([21 h 2 Y N+l (4/h 2 , h)Y N ^(4/h 2 , h) = AY N (A/h 2 , h) 2 ^ 0, and 
the result follows by induction. □ 

From Theorem 12.11 we know that Y n (t + h,h) — Y n (t, h) = hy n (t) + 0(h 2 ), and the 
zeroes of Y n (t, h) in t differ from those of y n {t) by 0(h 2 ), so (because the zeroes of y n 
are simple) we expect that Y n (t + h, h) should not vanish where Y n (t, h) does. Indeed 
this turns out to be the case. In the next section we shall see that the polynomials 
Yn are tau-functions for a discrete PII equation, and hence satisfy many other bilinear 
identities in addition to ( 1241) . To prove the following result we make use of one such 
identity, namely fl50l) below. 

Theorem 3.1 For any n = 1,2,..., if t is a zero of Y n , so that Y n (t ,h) = 0, then 
Y n (t ±h,h)^0 and Y n+1 (t Q , h)^0. 

Proof: It is true for n = 1 by inspection. Suppose the result is true for n = 
1, 2, . . . , iV — 1. Now assume that both Y/v(^o, h) = and Y^ito + h,h) = 0. Using 
the identity fl50l) with r n = Y n (t, h), and setting n — N — 1, t — to, we have 

Y N (t + h,h)Y N _ 2 (t ,h) - Y N (to,h)Y N _ 2 (t + h,h) 

= {2N-l)hY N ^{to + h,h)Y N ^{t ,h). 

By the assumption, the left hand side vanishes, and 2iV — 1 ^ for iVeZ which means 
that YAr_i(to, h) or YAr_i(t + h,h) = 0; but this contradicts the inductive hypothesis, 
so Y N (t Q + h, h) 7^ whenever Y N (t , h) = (and similarly for Y N (t — h, h)). Now if 
>at(£o, h) = 0, then from ( |24l) we see that 



h = 0. 




for all n G N. 



/T%-i(to, h)Y N+1 {t , h) = {h 2 to - 4)Fjv(t + h, h)Y N {t - h, h) 
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If we now assume that also Y/v+i(to> h) = 0, then the left hand side of the above vanishes, 
and since h 2 t — 4 ^ by Lemma [3TT1 then Y^(t + h, h) — or Y/v(^o — h,h) = 0, which 
is a contradiction. Therefore the result holds for n = N, hence for all n by induction. 
□ 

Clarkson and Mansfield have numerically studied the location of zeroes of y n (t) for 
low values of n [9]. They found that for a given n they occupy approximate triangular 
arrays in the complex plane and that the zeroes of y n {t) interlace in a certain sense 
with those of y n+ i{t). Since the zeroes of Y n (t, h) are the same as those of y n {t) up to 
0(h 2 ), the same picture holds for sufficiently small h, and our numerical observations 
show that the same qualitative behaviour persists for values of h up to at least order of 
unity. 



4. Alternate discrete Painleve II 



The second Painleve equation (Pn) is 

q = 2q 3 + tq + a, (27) 

where a is a constant. This second order differential equation is equivalent to a first 
order system introduced by Okamoto [14J, namely 



q = p-q* - 2 



,2 t 

(28) 

p = 2qp + £, l = a+\. V ' 

Rational solutions of ( |27l) arise when a = n G Z. However, in this and subsequent 
sections n need not be an integer unless stated explicitly otherwise, and we use the 
parameters n and a interchangeably. We also find it convenient to use the shifted 
parameter t = n + 1/2 as above, which fixes a point in the Ai root space. This is 
the parameter used by Okamoto to label solutions of Pu and the corresponding tau- 
functions, but we stick with the label n to maintain contact with the preceding results. 
Pn may be put into bilinear form by setting 

Then (12 71) with a = n is equivalent to the bilinear equations 

D\T n -T n _ x = F(t) T n T n _i, 



(D* - tD t + n)r n ■ r n _i = 3F(t) D t r n ■ r, 



(30) 



n— 1 ; 



where D t is the usual Hirota operator, and the function F(t) is arbitrary. By a gauge 
transformation on the tau-functions, that is r n — > G(t)r n for all n, the function F can 
be set to zero in (!30~!) without loss of generality. 

The Backhand transformation linking a solution of (1271) for a = n with that for 
a = n + 1 is 

q n +i = -q n , p n = q n + ql + t/2, (31) 

Pn 
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and the two solutions are related by 

Qn + ql = -q n +i + ql+i (32) 

(which can be seen as being inherited from the Miura transformation for KdV). From 
( pOj) and (13T1) it may be proved that, with the same choice of gauge that fixes F = in 
( I30|) . any three adjacent tau-functions r n , r n±1 satisfy the bilinear relations 

Ar„+i-r n _i = 2lr 2 n (33) 

and 

Tn+lTn-! = tr% - 2D\r n ■ r n . (34) 

With this choice of gauge, the equation ( 1341) is the bilinear form of the Toda lattice. 
As we have already mentioned (cf. equation ([3]) in the Introduction), the latter relation 
is the defining recurrence for the Yablonksii-Vorob'ev polynomials; with r n = y n they 
satisfy the other bilinear equations (133]) and fl30l) for F — 0, and hence determine the 
rational solutions of Pjj. 

Several discretisations of Pu have been studied in recent years including a discrete 
bilinear form [18] which tends to (1301) in the continuum limit. The discrete bilinear 
form allowed the construction of rational solutions to this discrete Pu, given in terms of 
ratios of determinants of Jacobi-Trudi type [19] , while solutions in terms of discrete Airy 
functions were found in [20, El] . A hierarchy of higher order analogues of this discrete 
Pu (dPu) have been given by Cresswell and Joshi, based on a Lax pair. However, due to 
the non-uniqueness of the discretisation process, there are many other dPII equations 
with analogous properties. For example, four different dPu equations are mentioned in 
[23], three of which are g-type (q-Pjj). Discrete Airy function solutions of one such q-Pu 
equation are found in [21], while another one is considered with its ultra-discretisation 
in [25], along with yet one more distinct q-Pu- In [26] a detailed study of a q-Pu in 
Sakai's classification [27] has been performed. Another alternative discretisation, which 
is relevant here, is the so called alternate dPu (alt-dPjj) which was studied in [12], and 
appeared in connection with the original dPjj in [28J. 

Here we consider the alt-dPn equation in the form of the second order difference 
equation 

%m-2 4r(m-l)-2 1 h 3 , 

^—T+ 2 ' = -9 + - + T (m-£)-2 (35) 

99 + 1 99 + 1 g 2 

with the notation g = g n (t), ~g n = g n (t + h,h), g = g n (t — h,h) 0. This is referred 
to as a discrete Pjj equation due to the fact that (1271) arises from it by setting 
g(t, h) = —l+hq(t, h) and taking the continuum limit h — > (with parameter a = £—1/2 
as usual). If we set g n (t, h) = x m with z m = h 3 m/2 — 2 and fi = —h 3 £/2 then this is 
equation (1.3) in [T2] (except that we have m in place of n). 

| The introduction of m in place of t is not essential, but it suggests that m and £ might be considered 
on the same footing, so that wg^ m = g n {t) should provide particular solutions of a suitable partial 
difference equation. See [T2] for a connection with the lattice modified BSQ equation. 
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The alt-dPjj equation (1351) has rational solutions for integer values of the parameter 
n = i — 1/2, provided by suitable ratios of the discrete Yablonskii-Vorob'ev polynomials 
Y n (t,h). In order to prove this, we must derive appropriate discrete analogues of the 
bilinear identities (!30|) . (|33|) and (|34|) . We start by presenting (!35j) in the form of a 
system, which makes it more manageable. 

Lemma 4.1 The alt-dPjj equation / T53]) for £ = n + ^,m = t/his equivalent to the 
first order system 

g n = (-2 + h 2 Pn )/(4: + 2g n -h 2 (g nPn + t)), 
Pn = (Pn-(hg n )/g 2 n . 

Proof: By rearranging the first equation in (1361) we have 



h 2 v ~L + g n g n 

and by shifting t — > t + h the latter provides a formula for p n in terms of ~g n and g n , and 
then substituting this into the left hand side of the second equation (I36I) . with p n given 
by (1371) on the right hand side, gives a relation between g n , g n and g n . After rearranging, 
shifting £ — > t — h, and setting t = mh, this relation is precisely (1351) . Conversely, given 
the alt-dPII equation (1351) . we can define p n by the formula fl37|) . which is equivalent 
to the first equation in ( 1361) : the second equation for p n then follows immediately from 
alt-dPIL □ 

The Backlund transformation for the alt-dPji equation (135]) is most easily derived 
starting from the analogue of (|32|) . and the following results are easily verified by direct 
calculation. 

Lemma 4.2 If g n +i is given by 

9n+i = 77-7^ > (38) 

Eh + g n p n 

in terms of g n and p n satisfying IfSh}) . then the identity 
1 1 

— + 9n = + g n +i (39) 

9n 9n+l 

holds. 

Corollary 4.1 The quantity g n+ i, defined by (E5|) with p n given by satisfies 
with I — > I + 1 . Equivalently, the equation / Tggj) and the relation 



I f (4-h 2 t)g n+1 n \ 

II \ 1 + 9n+l9n+l J 

together constitute a Backlund transformation for the alt-dPjj system IfSh}) . 

Corollary 4.2 The Backlund transformation for the alt-dPjj system, given by the 
formulae and has the following consequences: 

9n+lP„ = 9nPn] (41) 
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I ih 

— -0rH-l = — ! (42) 
9n Pn 

1 ih ( 1 \ 

£n+l#™ — = Sn + — J • (43) 



9n+l9 n Pn \ 9 

Remark. With g n — — 1 + hq n , the identity (132]) arises as the continuum limit of (139]) . 
and the formula ( 131"]) arises from ( 1381) . as /i — > 0. Similarly, the system ( 1281) is the 
continuum limit of (1361) . Nijhoff e£ al. derived equivalent formulae of Miura/Schlesinger 
type for the Backhand transformation of alt-dPu by making use of a variable y n (see 
equation (5.1) in [12]), which (modulo rescaling and replacing n by m) is analogous to 
p n defined by fl37|) . 

We now describe the tau-functions for the alt-dPu equation, which satisfy analogues 

of am. 



Proposition 4.1 Up to a choice of gauge, every solution of If35*\) is specified by a pair 
of tau-functions r n (t, h), r n ^i(t, h) via the formula 

*.(«.*) =- T -'L "'Ti't w 

r n _i(*,/i)r n (t - h,h) 
where the tau-functions satisfy the bilinear equation 

TnZn-l + ZrJn-l = 2r n T n _i (45) 

and the quadrilinear ( degree four) equation 

(4 - mh^r^T^Tn^^ - T n _{r n ) + (4 - (m - l)/i 3 )r n _ir„(7:„_i^ - Z n _{r n ) ^ 

+8(ZI-I r n - T n-lln) ~ ^nh 3 '' T n -lZn T n-lT n = 0, 

m = t/h, n = a = £ — 1/2 and r n = r„,(t, /i),r n = r n (t + h, h),T n = r n (t — h, h), 

etc. 

Proof: Upon substituting the tau-function expression (144")) into (|35l) and clearing 
denominators, a relation of degree eight results, which can be simplified somewhat 
by rewriting it in terms of the symmetric/antisymmetric quadratic quantities A± = 
r„_iT„ ± r„_,T n and A, = r n _iT ± r t„. In general, for any choice of tau-functions 
the bilinear equation 

TfiZn— 1 ~^~Zn T n-l = T n T n _i 

holds, for some function F = F(t, /i), but by applying a gauge transformation r n — > G r n , 
r n _! — > with GG/G 2 = F the function F can be removed to yield the bilinear 

equation (1451) . With this choice of gauge, the remaining terms in the degree eight relation 
factorise to yield the quadrilinear equation (|46p . and conversely if these two tau-function 
equations hold then g n given by (jUJ) is a solution of (j551) for £ = n + 1/2 = a + 1/2. □ 



Remark. The existence of a quadrilinear relation between a pair of tau-functions is 
mentioned in section 4 of [T2|, where a third tau-function is introduced to obtain purely 
bilinear relations (cf. Theorem 14. II below) . 
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It is easy to see that (I45p tends to the first of (15U|) in the continuum limit (with 
the gauge chosen so that F — 0). Although the second relation (T461, between the two 
tau-functions is of overall degree four, it still produces the second bilinear differential 
equation (1301 in the continuum limit (provided that the first one also holds). In order to 
work with purely bilinear equations in the discrete case, we must consider three adjacent 
tau-functions r n , T n ±\. 



Theorem 4.1 Up to a choice of gauge, every solution of the alt-dPjj system ( fjpj) is 
specified by three tau-functions r n _i(i), T n (t), r n+ i(t), with g n given by ( fj^p and 

where the tau-functions satisfy ( fj^ ) as we// as 

^n+lZ„ + Zn+l T n = ^T n+1 T n (48) 

and 

/i 2 r n+1 r n _! = (h 2 t - 4)r nZn + 4r 2 . (49) 
With this choice of normalisation the identities 

r n+l r n-l — T n+l r n-l = 2£hT n T n (50) 

and 

Tn+lZ^ ~ -Ln+^n-l = 4lh,T* (51) 

also hold. For this choice of gauge, these purely bilinear relations are compatible with 
the Bdcklund transformation ( fff£j) for alt-dPn, in the sense that r n and r n+ i satisfy ftjty ) 
for n — > n + 1, and 

/. t \ _ r n (t - h,h)r n+1 (t,h) 

9n+l(t, tl) — T—rz 77 —rr 

T n {t,h)T n+ i(t - h,h) 
satisfies (E^j with £—>£ + l. 

Proof: If a solution g n , p n of ( 1361) is given by the expressions ( )44l and (H7J) respectively, 
and the gauge is fixed by f|45|) . then the latter implies that 



1 r 2 1 
+ — = -2^^ = + g n+1 

9n 7~nT-n 9n+l 

by Lemma [4.21 where g n +x given by (1381) is a solution of (1351) with £ —>■ £ + 1. The first 
equality above implies fT49|) . while the relation fj4Tl) implies that (? n+ i is given in terms 
of tau-functions by the formula fj44l) with n — > n + 1, and hence fT48l) follows from the 
second equality above. The bilinear identities ( 1501) and ( 1511 then hold as a consequence 
of the relations (j42|) and (jUJ) respectively. By Proposition 14.11 the given choice of gauge 
implies that the pair r n _i, r n satisfy ( 1461) . and the fact that (? n+ i is a solution of alt-dPII 
with the parameter shifted implies that the pair r n , r n+1 also satisfy this quadrilinear 
equation with n — ► n + 1. 
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Conversely, suppose that g n is defined in terms of tau- functions by (pHj) . g n+1 is 
defined by the same relation for n — > n + 1, and p n is defined by (1471) . where the tau- 
functions satisfy the three relations ( )45j) . ( l4~8]) and ( 149]) . The identities ( 1371) and ( 139]) 
follow immediately. Furthermore, from (145]) and (148]) it is clear that 

r n-l / r n+l , Hn+1 \ _ r n+l / r n-l Z n -1 



which implies that 

T~n+lT~n-l ~ Tn+lTn-1 _ r n+lZ n -l — Zn+1 T «-1 

Therefore (T n+1 r n _! — r n+ iT„_i)/ (2hr n T n ) is independent of t, and if we denote this by 
£, then we have the bilinear equation floU]) . which implies that f]4"2~]) also holds. Solving 
(142]) for g n+1 gives an expression in terms of p n and g n , which in turn means that g n+1 
can be written in terms of g n and ~g n using (13 7p . By shifting t t — h, this gives a 
formula for g n+1 in terms of g and g n , and then substituting for g n+i and g n+ i in (139]) 
yields the alt-dPn equation ( 135]) . It then follows that (? n ,Pn satisfy the system ( |36]) . □ 

Theorem 4.2 For parameter t = n + | roi/i n G Z i/ie alt-dPu equation [3b^l has 
rational solutions given in terms of the discrete Yablonskii-Vorob'ev polynomials by 

Y n ^(t-h,h)Y n (t, h) 
9n ~ Y n ^{t,h)Y n {t-KhY 
where the polynomials are extended to negative indices n by setting YL n = Y n -\ for 
jigN. As well as the defining recurrence [24\ ), the relations (4^ , [Jo] ), / f50j) and 137]) 
are satisfied by r n (t, h) = Y n (t, h) for all n G Z. 

Proof: When i = 1/2, the equation (]35j) has the trivial constant solution go(t, h) = —1, 
which can be obtained by setting r = Y = 1 = F_i = r_i in ( |44"j) . and from (|37j) we 
have po = £/2 which gives ri = Y\ — t by (l47j) . It is easy to verify that each of the 
bilinear equations ( 14"5]) . ( l4~8j) and ( H9]) is satisfied by these tau-f unctions. By applying 
the Backlund transformation (138]) repeatedly (both forwards and backwards) a doubly 
infinite sequence of rational solutions { g n } ne z is obtained. Then by Theorem 14.1] 
since the Backlund transformation is compatible with the choice of gauge, it follows 
by induction that the corresponding tau- functions satisfy the identities (]24l) . (]4"5]) . f]46]) . 
( 150]) and ( |5T]) for all n G Z. Since the Yablonskii-Vorob'ev polynomials are defined by 
(149]) with r = 1, r x = t, it follows that this particular sequence of tau-functions is 
given by r n (t, h) = Y n (t, h) for all nGff. The fact that this relation can be consistently 
extended to negative n follows from the observation that all of the tau-function identities 
in Proposition 14.11 and Theorem 14. II are invariant under n — > —n — 1, £ — > —I. □ 

Remark. The simplest rational solutions of alt-dPu (corresponding to n — 0,±1) 
are described in section 6 of the paper [12] by Nijhoff et ai, where it is indicated how 
the above sequence of rational solutions can be generated recursively via the Backlund 
transformation, but no closed form for these rational solutions is given in that work. 
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The fact that the alt-dPu equation can be derived from a sequence of Backlund 
transformations applied to solutions of Pjjj provides a relation between the discrete 
Yablonskii-Vorob'ev polynomials and the Umemura polynomials for Pni- The third 
Painleve equation, Pm, is given by 

,'\2 „„/ 1 X 



[W 



w 



w 



— + -(aw 2 + (3) + -fw s + -, 



(52) 



W X X w 

with the prime ' denoting d/dx, and without loss of generality (by rescaling the 
independent variable x) if 7 5 7^ the latter two parameters can be fixed as 7 = —5 = 1. 
(Note that the parameter a in (152]) should not be confused with the parameter a in 
P/j.) Backlund transformations for Pjjj (cf. section 2 in [T2]) can be used to relate 
three adjacent solutions w n = w(x; a,/3) and io n ±i — w(x; a ± 2,(3 ± 2) and if we set 
g„ = —l/w n then this contiguity relation can be written in the form of the alt-dPjj 
equation 



Zn-1 



9n+l9r. 



9n9n-l 



where 



-T + 7: \9n + — ) -n I' 0. 

12V 9r, 



(53) 



z n = (a + (3 + 2)/4, Az n := z n+1 -z n = 1, [i = (J3-a- 2) /A. (54) 

(To compare with equation (2.5) in [T2], set g n = ix n , x = t in the above.) For a certain 
set of parameter values, Pm has rational solutions which are described by the following 
result. 

Theorem 4.3 (Kajiwara & Masuda [73]/ ) For parameters 

a = 2n + 2v - 1, /? = 2n - 2v + 1, neZ (55) 

and 7 = —5 = 1, the third Painleve equation / f5^j) /ias rational solutions w(x; a, (3) = w n 
given by 

V n (x,u- l)V n _i(x,u) 



T> n (x, v)V n _ x (x,v- 1) 
where the polynomial T> n is given by a determinant of Jacobi-Trudi type, 



(56) 



V n (x, v) 



Pn Pn+1 
Pn-2 Pn-1 



■ ■ Pln-\ 

■ ■ P2n-3 



P- n+ 2 P-n+3 •■■ Pi 

with p k = Pk(x, v) defined by the generating function 



^2p k (x,v)\ k = (1 + A)^exp(xA) 



/c=0 



and p k = for k < . 
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Remarks. The polynomials pk{{x,u) are essentially just associated Laguerre 
polynomials, and, for each n, T> n (x, v) is a Schur polynomial with restricted arguments, 
corresponding to the partition (n, n — 1, . . . , 2, 1). The result stated above is an adapted 
form of Theorem 1 in [13], and describes one family of rational solutions of Pi if, for 
a complete description of all rational solutions of Pm for 7$ 7^ 0, see [29J. The 
polynomials T> n (x, v) (after some scaling) are known as the Umemura polynomials 
for Pm. Further properties of scaled Umemura polynomials for Pm are detailed in 
[30J, including the remarkable patterns formed by the roots, and differential/difference 
equations; analogous polynomials corresponding to the special cases when 7<5 = are 
also treated there. 

Theorem 4.4 The discrete Yablonskii-Vorob'ev polynomials are given in terms of 
determinants of Jacobi-Trudi type by the formula 

Yn(t,h) = Cn h^^V n ^l-±y (57) 

where 

c n = (2n-l)!!(2n-3)!!...3!!l!! (58) 

for nGN. 

The proof of the preceding theorem makes use of some results in the next section, 
and is relegated to the appendix. However, it is clear that if we rearrange the formula 
( |57|) then we can rewrite T> n (x, v) in the form of a Hankel determinant. 

Corollary 4.3 The Umemura polynomials for Pm, given in scaled form by TJ n (x, v), 
are proportional to the Hankel determinants u n as in ( 173]) with h = (x/4)~ 1//3 and 
t = ( x /4:)- 1/3 (x + v) . 

Remark. It is known that the Painleve differential equations form a coalescence cascade 
from Pyi down to Pi (see J3T]). In [13] it is shown that the coalescence limit from Pm 
to Pa produces the Yablonskii-Vorob'ev polynomials y n {t) as a limit of the Umemura 
polynomials, but this arises in a different way compared with the limit h — > considered 
above. More precisely, with the scaling used here, the coalescence from (152]) to (127]) arises 
when the independent variable x and parameter v scale asx=^ + ^-,z/ = | — -5, with 
e — > 0. In this limit, up to scaling the polynomials V n (t/e + 4/e 3 , 1/2 — 4/e 3 ) produce 
y n (t) at leading order in e. 



5. Symplectic properties and discrete Pxxxiv 

Okamoto [H] showed that Pn can be written as the system (I28p which is in Hamiltonian 
form, i.e. 

dH dH 
?=-5-> P = - — , (59) 
op oq 
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with 

H = ^-(q 2 + t -)p- Iq. (60) 
Eliminating p gives ( 1271) . which is Pn, whilst eliminating q gives 

p=^ + 2p 2 -tp-— (61) 

which is known as Pxxxiv (see [31], chapter XIV). 

This representation has been used in [32] to study further properties of the 
Yablonskii-Vorob'ev polnomials. Although the alt-dPn equation ( l35l) . being a non- 
autonomous difference equation, does not have a Hamiltonian form, many of the results 
proved there have their counterparts in the discrete case. For ease of comparison, we 
briefly recall some known results on Pjj. In terms of the canonical coordinates (q n ,Pn), 
the Backlund transformation for Pjj can be written (together with its inverse) as 

q n+ i = -q n - Pn+i = ~Pn +t + 2(q n + ^ 

g"-l = -<ln- 2£-p n +V Pn-1 = -Pn + t + 2q 2 n . 

It is straightforward to check that dq n+ i A dp n+ \ + dH n+1 A dt = dq n A dp n + dH n A dt, so 
the transformation n — > n + 1 is a canonical transformation on the extended phase space 
with coordinates (q n ,p n ,t). The generating function for this canonical transformation 



(62) 



is 



so that 



2 

F{qn, Qn+i) = I log(g„+i + q n ) + -ql +l + tq n+1 (63) 



Pn — ~ T\ , Pn+1 



oq n oq n+1 
The Backlund transformation formulae ( 1621) imply that any sequence of solutions q n of 
Pn (labelled by n = £ — 1/2) satisfies 
f £ — 1 

+ + 2q 2 n + t = 0, (64) 



Qn+1 + Qn Qn + q n -l 

which is a discrete form of Pj, whilst p n satisfies 

(Pn+1 ~ Pn-lYp n ~ ^{Pn+1 + 2 Pn + Pn ^ - 2t)p\ + 4£ 4 = 0. (65) 

We will now show how the above results carry over into the discrete case. 

Proposition 5.1 In terms of the variables q n = (g n + l)/h and p n , the Backlund 
transformation for the alt-dPjj equation K35\) ( corresponding to n — > n + 1 ) can be 
written together with its inverse (corresponding to n — > n — 1) as 

q ' n+1 h(q nPn +e)-p n ' 

Pn+i = (q n Pn + £f Q| - - £(* - 2p n )(q nPn + £)-p n + t; 

_ __ I l-i+hq2-2ql+q n p n -q n t+l-l)+h 2 ql{-2p n +t)+h- i qip n ^ > 

Qn-l qn T - Pn+ t+2qi-hq n {2ql-Zp n +2t)+Wql(-?, Pn +t)+hZqlp n ' 

Pn-i = ~Pn + t + 2ql + h(2p n -t)q n - h 2 p n q\. 
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This Backlund transformation is a symplectic map in the phase space with coordinates 
(q n ,Pn), with generating function 

*.+>) = %^ + to - U hH ~ 4 + ehS ) ~ 

+ t \og(hq n q n+1 - q n+1 - q n ) (67) 

such that dJ- = p n +idq n +i — p n Aq n ; in other words the canonical Poisson bracket 
{qn,Pn} = 1 is preserved for n — ► n+1. The equations KUh}) and the generating function 
reduce to those of the Backlund transformation for Pjj in the continuum limit h — > 0. 

Proof: The formulae fl66l) follow immediately from Lemma 14.21 and Corollary 14.11 To 
verify that the map (q n ,p n ) i— > (q n+ ±,p n+ i) is symplectic, it is sufficient to calculate 
directly that its Jacobian determinant is equal to 1. This also follows directly from the 
closure of the exact one- form dJ- ' upon checking that p n = — p n +i = a 9:F with the 
generating function T as in the formula (1671) . It is straightforward to verify that the 
relations fl66l) have the correct continuum limit given by fl62|) . and that flBTl) also yields 
(1631) when h -»• 0. □ 



Remark. The symplectic structure for the alt-dPjj equation can be derived from 
Okamoto's Hamiltonian formulation for Pru, since the canonical coordinates (q n ,Pn) 
above are related to those for Pj n by a shift and rescaling (cf. section 2.3 in [30] for 
instance). 



Corollary 5.1 For the alt-dPjj equation (fff5j). the analogue of (64\ ) is the equat 



ion 



' C- 1 * + T ^L+ t - M = . (68) 



q n +i + q n ~ hq n+1 q n q n + g n _i - hq n q n _i 1 - /ig ; 
In terms of g n this is the second order difference equation 

' +T <^ + iL( fc + iW<«=o, (69) 



SWflSVi - 1 0n0n-l ~ 1 /i 3 \ h 3 h 

which is another form of the alt-dPu equation. The conjugate momentum p n satisfies 
a third order recurrence relation in n, namely 

f p n _ fePn[(2n-l) 2 +2p^_ 1 (p n -p»_ 2 )] _ h 2 p„t _| ftV(2n-l)i 
Pn+l\ 2n+ i 4p n „i(4n 2 -l) 4(2n+l) 8(2n+l) > 

I r Pn-1 _ fep, 1 -l[(2n + l) 2 +2p 2 1 (pn-l-Pn + l)] _i_ ft 2 Pn-lt j_ fe 3 p n - 1 (2n+l) -[ 

T-^n-21. 2n-l 4p„(4n 2 -l) " l " 4(2n-l) ^ 8(2n-l) J 

, fe 2 PnPn-i , fe(4n 2 -l) _ (2n+l)[fe 3 (2n-l)-2fe 2 f+8] (7Q\ 
~+~ 2(4n 2 -l) + 8p„p„_i 16p„ V ; 

(2n-l)[fe 3 (2n+l)+2fe 2 f-8] . r 8-2fe 2 t+(4n 2 -3)fe 3 1 _ (6ra+5)fep„ 

16p n _! +PnPn-l[ 4(4n 2 -l) J 4(2n+l) 

- (6w 4 ~ 5 ^- 1 + A[ 3 2t + I6/1 - AhH 2 - 4/i 3 t + (An 2 - l)h A ] = 0. 

Proof: Upon solving the first of ( l66l) for p n and substituting in the last one, (JB"8l) results. 
One can eliminate the quadratic terms in q n from the second and fourth of (1661) to get 



2n + 1 2(p n+ i - p n _i)p n + (t - 2p n )(2n + l)h 
qn ~ Ap n 2(2n + l)(h 2 p n -2) [ ' 

After substituting this into the first of (|66|) with n — > n — 1, we obtain (I7D1 . 



Rational solutions of discrete Toda and alt-dPjj 



18 



Remarks. The fact that the alt-dPjj equation is self-dual, in the sense that the 
superposition formula fl69l) for its Backund transformation is (up to rescaling and 
reversing the roles of the dependent variable and the Backlund parameter) the equation 
itself, was first noted in [12]. Equation (37) in [32] is the continuum limit of equation 
( 1701) . The latter relation allows one to obtain p n +i uniquely given p n ,Pn-i,Pn-2- 
Remembering that p n = Tn - 1 ^ n+1 ; W e see that if Tj(t, h) for j = n — 3, n — 2, . . . , n + 1 
are given for a particular value of t , then r n+2 can be evaluated for this same value 
of t. More precisely, (1701) is equivalent to a recurrence relation for the r n that involves 
no shifts in t and is linear in r ra+2 and r n _ 3 ; in particular this relation, which is of fifth 
order in n, is satisfied by the discrete Yablonskii-Vorob'ev polynomials Y n for n G Z. 
The latter recurrence for the tau-functions, which is omitted here, tends to equation 
(38) in [32J in the continuum limit. It is also possible to write q n in terms of unshifted 
T n , by substituting the right hand side of (HTj) for p n in ( I7T1) . 

We now consider the map in the (q n ,Pn) phase space corresponding to shifting in t 
rather than n. 

Proposition 5.2 In the phase space with coordinates (qn,Pn)> the shift t — > t + h 
corresponding to the alt-dPn system [36]) is given by 

— _ -2qn+h(t—2pn)+h 2 p n qn 
In — -2-2hq n +h 2 {t-p n )+h :i q nPn i 

Pn = (fc ^_ 2) 2 ( - 2 - 2h 1n + h 2 (t- p n ) + h 3 q n p n ^j X (72) 

x ( - 2p n - 2h(£ + q n p n ) + h 2 p n (t - p n ) + h 3 p n (£ + q n p n 
This is a symplectic map with the generating function 
S = ^(A-hH)\og(2-h 2 p n ) + lhg^ + l(p n -p^ 



+ + Ap n p n - Itantr 1 ( > 2 ^ 4 ^ 

such that dS = q n dp n — q n < 



(73) 



(74) 



Proof: The equations (1721) just correspond to rewriting ( |36l) in terms of q n = (g n + l)/h, 
q n = (g n + l)/h. It is extremely easy to check directly from fl36|) that the symplectic 
form u n = dq n A dp n = \dg n A dp n is preserved by the shift in t. To find the generating 
function, it is convenient to write g n ,g n in terms of p n ,P n i giving 

H = 9n = ~ (£&) ~ ifc Q =F + A VnPn ] 

(One has to take the upper choice of sign in each case to get the correct continuum 
limit.) Having found an 5* for which the relations (1741) hold, S = (S + p n — p n )/h 
provides the generating function in ( J73l) . □ 

Corollary 5.2 The discrete Pxxxiv equation associated with the alt-dPjj equation If3h]) 
can be written either as 

± v^ 2 + ApnP n ± JPh 2 + Ap n p n = 2P ;^ ~ hH) , (75) 

a n p n 
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or with the square root signs removed as 

(>/> 2 + 4 Pn p n ) + ±VnV_ n ) = ( 2 fl^j ~ " 2MVn + 2 ■ 

Proof: The second order recurrence relation for p n is obtained by downshifting the 
second of (|74p and equating it to minus the first. □ 

Remark. The 3-point correspondence (175]) is equivalent to an analogous equation 
for the variable y n , that is equation (5.3) in [12], and the equation (175"]) has the 
same structure as certain discrete Ermakov-Pinney equations constructed in [33]. The 
presence of this structure is due to the connection with discrete Schrodinger equations 
(for which, see the proof of Proposition 16. II below) . 



6. Lax pair for alt-dPjj 



The continuum PI I is equivalent to the pair of bilinear equations (I30j) . A Lax pair for 
Pn is given by the linear problem 



where 



and 



m 

~dt 



B 



_ZZ 

2 

iTn-1 
2r n 



dr] 



2 

2 



A*, 



A 



-v 2 -l + 

iVTn-l _|_ 

T n r dt 



Tn-lTn + 1 

2T 1 

Tn 



T„(t) ' dt \ T n 



Tn-lTn + 1 

2r„2 



Consistency of (jTbj) requires 
dB _ dA 
dr] dt 

leading to the two conditions 



+ [B,A] = 0, 



<t>± + Vfa 



0. 



(76) 



(77) 



(7f 



(79) 



10) 



where we have set 



T 27V, ' 



V 



t 

2 



The choice of gauge V = t/2 — \ 



2-C- logr n gives precisely the equation 



'dt 



and 



with this normalisation for the tau-functions the conditions (180{) give the first equation 
in (I3"U|) f or F = 0, together with the same equation for n — > n+1. The bilinear equation 
(133]) is a consequence, and the second equation in (1301) then follows. It is well known 
that (179]) is an isomonodromy condition: the monodromy of the solutions of the second 
linear equation ( 1761) in the complex 77 plane is independent of t if and only if Pn holds. 

In the discrete case the situation is completely analogous, based on a linear problem 
that comes from the first part of the discrete Toda Lax pair ([5]). 
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Proposition 6.1 The linear problem 

* = BV, = Am (81) 



where 

and 
A = 



b=\ _L T ') = (i*A ' (82) 



/; :; \ A 2 



! + if 4 -2(0>_ + l) + A 2 H ' 1 



constitutes a Lax pair for the alt-dPu equation /[35)) , which is equivalent to the 
consistency condition 

\d\B + BA-AB = (83) 

for ( [Hp . In £/ie /zmzi /i — > 0, this linear system reduces to the Lax pair [76] ) for Pjj. 

Proof: The consistency condition (I8"3"|) implies the two relations 

0± + ^± + ± = O, (84) 
which take the form of discrete Schrodinger equations with 
v _ (hH-A) _ 2r n 2 (/A-4) 

2(1 - 0+0_) AtI - /i 2 r n+ ir n _i ' 

If we fix the gauge so that V = — 2r 2 /(r n r n ), then we get precisely the relation f|49|) . 
and upon substituting the latter expression for V into each of the equations ( |84|) in turn 
we find that fl4"5"j) and (|4"8"1) hold. Then by Theorem 14.11 g n = — T n _iT n /(r n _ir n ) satisfies 
the alt-dPu equation. For the continuum limit one should take A = , which gives 

A = ^A + O(h ),B = l + hB + 0(h 2 ), (85) 

so that the condition f!79|) arises from (1831) as /i — > 0. □ 



Remark. In [12] a different 2x2 Lax pair is presented for the alt-dPu equation, by 
reduction from the modified Boussinesq lattice. However, we have not found a direct 
relationship between these two Lax pairs. 



7. Concluding remarks 

There are many ways to construct a discretisation of a given integrable differential 
equation, depending on which particular properties (e.g. Lax pair, explicit solutions, 
Poisson structure, Hirota bilinear form,. . .) one most wishes to preserve. (For a 
thorough account of the Hamiltonian approach, see [34] .) Due to the non-uniqueness of 
discretisation, it is not always clear what is the "best" discrete analogue of a continuous 
system. The derivation of the alt-dPu equation presented here was initially motivated 
by the construction of exact rational solutions, but it has turned out that analogues 
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of all the other structures associated with the second Painleve equation arise naturally 
here as well. The fact that this discretisation scheme for Pjj, based on the discrete 
Toda lattice, turned out to be connected with the superposition formula for Pm was 
completely unexpected by us, but led to different expressions for the polynomial tau- 
functions in terms of Jacobi-Trudi determinants. In [T2] other solutions of alt-dPn were 
constructed in terms of Casorati determinants of discrete Airy functions. In future we 
should like to analyse other solutions of the equation fl35|) . We constructed the rational 
solutions of this equation from polynomial tau-functions given by Hankel determinants, 
but recent results for the continuous case [7] lead us to expect that all tau-functions 
should have a similar structure. 
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Appendix 

Here we present the proof of Theorem 14.41 The result essentially follows from the fact 
that the Backhand transformation for Pm generates a sequence of rational solutions w n 
given by ( 1561) . which simultaneously provide rational solutions of the alt-dPu equation 
f l53|) by setting g n = —l/w n with z n = n + 1/2 for n G Z and \i = —v. On the other 
hand, Theorem 14.21 and Corollary 15.11 together imply that the alt-dPn equation in the 
form ( 1691) has rational solutions, given by suitable ratios of discrete Yablonskii-Vorob'ev 
polynomials, when i = n + 1/2 with nsZ. Upon comparing (|53l) with (|69j) we see that 
these two sets of rational solutions coincide if we identify x = 4//i 3 , v = t/h — A/h 3 , and 
then it follows from Theorem 14.31 and Theorem 14.21 that 



for all n G Z, where Z n (t,h) denotes the right hand side of fl57|) for n > 0, and the 
identity extends to negative n upon setting Z_ n = Z n _\. Then it is necessary to show 
that Y n (t, h) = Z n (t, h) for all n. It is sufficient to consider n G N as in Theorem I4.4[ as 
the extension to negative n is trivial. 

For n = the result Y = 1 = Z is obvious, so to use induction we assume that 
Yat_i = Zn_i and then from ( 1 A. II) with n = N it holds that 



Both sides of the relation (I A. 21) are ratios of polynomials in t (for Z^ this follows from 
the determinant formula for IV), and the rational functions on each side must have the 



Z n _x(t, h)Z n (t-h, h) = Y^jt, h)Y n (t-h, h) 
Z n _ x {t - h, h)Z n (t, h) Y n -i(t - h, h)Y n (t, h) 



(A.l) 



Z N (t - h, h) _ Y N (t - h, h) 
Z N (t,h) ~ Y N (t,h) 



(A.2) 
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same zeroes and poles. However, although Theorem 13.11 implies that the numerator and 
denominator on the right hand side have no common factors, we cannot immediately 
assert that the same is true on the left hand side without knowing the degree of T>n(x, v) 
in v. (The proof of Corollary 2 in [13] just gives the degree in x of V N , denoted F N there, 
as N(N + l)/2.) However, by Proposition 3 in [13J these Jacobi-Trudi determinants 
satisfy the recurrence 

(2n + l)P n+1 P n _! + x{V n V" n - (V' n ) 2 ) + V n V' n - (x + v)V\ = 0, (A.3) 

with V-i = 1 = V and ' denoting d/dx. In fact, this recurrence can be used to show 
that V N (x,v) is also of degree N(N + l)/2 in u, but we do not need this. Instead, 
setting v = in ( 1A.3I) leads to the expression T> n (x, 0) = c~ 1 x n ^ n+1 ^ 2 where c n is 
given in terms of double factorials by (1581) . Hence for the case at hand we have 
Z N (A/h 2 ,h) = c n h ni - n+i y 2 V n (A/h 3 ,0) = (4//i 2 )™(™ +1 )/ 2 = Y N (A/h 2 ,h) by Lemma O 
Using ljA.2l) it then follows by induction that Z N (A/h 2 + kh, h) = Y N (A/h 2 + kh, h) for 
all HZ, and since these two polynomials agree for infinitely many values of t they 
must be equal, as required. □ 
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